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Abstract 

We study the quark mass expansion of the magnetic moments of the nucleon in a 
chiral effective field theory including nucleons, pions and delta resonances as explicit 
degrees of freedom. We point out that the usual powercounting applied so far to 
this problem misses important quark mass structures generated via an intermediate 
isovector Ml nucleon-delta transition. We propose a modified powercounting and 
compare the resulting chiral extrapolation function to available (quenched) lattice 
data. The extrapolation is found to work surprisingly well, given that the lattice 
data result from rather large quark masses. Our calculation raises the hope that 
extrapolations of lattice data utilizing chiral effective field theory might be applicable 
over a wider range in quark masses than previously thought, and we discuss some open 
questions in this context. Furthermore, we observe that within the current lattice data 
uncertainties the extrapolations presented here are consistent with the Pade fit ansatz 
introduced by the Adelaide group a few years ago. 
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1 Introduction 



The computation of nucleoli properties in Lattice QCD is progressing with steadily increas- 
ing accuracy [[!]]. So far, these results are, however, limited to relatively large quark masses 
and the quenched approximation. The typical "light" quark masses manageable on the lat- 
tice up to now, are more than 10-20 times larger than the average u- and d-quark masses, 
m q ~ 8 MeV, determined at a renormalization scale around 1 GeV. This corresponds to 
pion masses well above 0.5 GeV. 

On the other hand, in the chiral limit m q — > 0, QCD at low energies is realized in the form 
of an effective field theory with spontaneously broken chiral symmetry, with massless pions 
as the primary active degrees of freedom. Nucleons are added to this theory as fermionic 
matter fields which can be treated non-relativistically due to their large mass. The coupling 
of the chiral Goldstone bosons (i.e. the pions) to these spin 1/2 matter fields (i.e. the 
nucleons) produces the so called "pion-cloud" of the nucleon, an important component of 
nucleon structure at low energy and momentum scales (for a survey and references see 0]). 

The interpolation between lattice QCD results and the chiral limit, passing through the 
physical pion mass point at ss 0.14 GeV, is presently a subject of intense studies and 
lively debate. Leinweber et al. || have initiated such considerations by requiring that any 
interpolation of this kind should be subject to the leading dependence on the pion mass, as 
dictated by chiral symmetry. 

In the present workQ we focus on the chiral aspects of nucleon magnetic moments. The 
authors of Ref.0 have introduced a parameterized form of the pion mass dependence of 
proton and neutron magnetic moments in order to interpolate between available (quenched) 
lattice QCD results and the "physical" values at the proper m n . Our approach uses in- 
stead systematic methods of chiral perturbation theory (ChPT) in the baryon sector 0. It 
will be demonstrated that such a scheme provides an interpolating function of m^, which 
successfully connects lattice data with actual physical magnetic moments, in such a way 
that the physics behind the chiral magnetic structure of the nucleon can be identified. Not 
surprisingly, virtual excitations of the delta isobar turn out to play a decisive role in this 
context. The strong N — ► A magnetic dipole transition induced by the photon field is 
a prominent feature in determining the non-linear dependence of the magnetic moments 
on the pion mass. Given that delta excitations are separated from the nucleon mass by 
a "small" scale of hardly more than twice m n , it is obviously not a good starting point 
to relegate the important physics of the A(1232) to higher order "counter terms", if one 
is interested in the m^-dependence of nucleon structure properties above the physical pion 
mass. Such effective chiral field theories with explicit pion, nucleon and delta degrees of 
freedom and a systematic power counting already exist in the literature (e.g. see Ref.|| and 
references therein). Here, however, we argue that even these existing schemes, such as the 
so-called "Small Scale Expansion" of Ref. |J, need to be modified if one wants to capture 
the important A(1232) induced m^-dependence already in a leading one-loop calculation. 

This paper is organized as follows. The chiral effective field theory framework, including 
the treatment of nucleon and delta degrees of freedom, will be summarized in section |2|. 

1 Some aspects of this work have already been reported in Rcfs.pL 0. 
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The actual calculation of the anomalous magnetic moments is described in section |3], with 
the results discussed in sections f|, ^] and ^|. The following two sections examine higher 
order corrections to be compared with the Pade-approximants suggested by the Adelaide 
group |J, and discuss some aspects of the quenched approximation as they may be relevant 
in extrapolations toward the presently existing lattice data. Undeniably, uncertainties still 
remain in such extrapolations to pion masses as large as five times the physical one. The 
results are nevertheless promising, and further progress in lattice QCD towards smaller 
quark masses in partially-quenched or even unquenched simulations can be expected to 
reduce the errors substantially. 



2 Effective Field Theory Input 
2.1 General Remarks 

The effective field theory that approximates QCD in the low-energy limit is Chiral Pertur- 
bation Theory (ChPT) J7|, [J. In the case of two light quark flavors its primary degrees of 
freedom are pions as Goldstone bosons of spontaneously broken chiral symmetry. ChPT 
represents a systematic expansion in terms of low-momentum and small quark mass scales, 
which allows for a systematic treatment of the additional explicit breaking of chiral symme- 
try responsible for non-zero pion mass. Baryons are introduced as fermionic matter fields, 
acting as sources which create or annihilate pion fields in accordance with chiral symmetry. 
In the non-relativistic versionP] of baryon ChPT considered in this work ("heavy baryon 
ChPT") 0, the nucleon mass M is treated as a "large" scale that persists in the chiral 
limit, allowing for a 1/M expansion in the theory. 

The A(1232) isobar is the lowest spin-3/2 excitation of the nucleon, reached via a strong 
magnetic dipole transition. It is therefore of key importance in considerations of the nu- 
cleon's magnetic structure. Since its mass differs from that of the nucleon by less than 0.3 
GeV, its incorporation as an explicit degree of freedom in baryon ChPT turns out to be 
mandatory in the present context. The basic techniques to do so are well established, at 
least in the context of non-relativistic baryon ChPT ||. However, we emphasize that a 
chiral effective field theory that couples both octet and decuplet baryons to the Goldstone 
boson dynamics in addition must specify how the small but always finite octet-decuplet 
mass-splitting — throughout this work denoted by the parameter A — is merged with the tra- 
ditional powers-of-g counting of ChPT spelled out in Ref.||. In that sense we state that 
a unique extension of spin 1/2 baryon ChPT to a chiral effective field theory with explicit 
spin 1/2 and spin 3/2 baryon degrees of freedom does not exist. An additional piece of 
information concerning the proper (ac) counting of A has to be provided, without guidance 
from chiral symmetry constraints^. Here we specifically address this issue by adopting the 

2 Some results of relativistic baryon ChPT are discussed in appendix [b|. 

3 For example, we note that a combined treatment of non-relativistic baryon ChPT and large iV c -counting 
rules can provide a systematic power counting for systems with both spin 1 jl and spin 3/2 degrees of freedom 
0- 
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(phenomenologically motivated) philosophy, spelled out in Ref.[ll], to count in the param 



eter e = (q, A) which, in addition to the (usual) powers-of-g, also keeps track of the finite 
mass-splitting A. All discussions in the subsequent chapters refering to leading order (LO), 
next-to- leading order (NLO), etc. Lagrangians are then understood as "powers-of-e" , i.e. 
0(e), 0(6 2 ),... 1. 

The following subsections give a brief summary of our basic input. We will work with 
2-flavor ChPT coupled to spin 1/2 and spin 3/2 matter fields throughout. We also note 
that strictly speaking all parameters mi, fi, ■ ■ ■ in the Lagrangians to be introduced in 
section |2| should carry an extra superscript index g®, m' ... to distinguish them from any 
physical quantities which might carry the same label. In most cases we suppress this extra 
index to obtain simpler formulae, but it is understood that there is a difference, even if it 
sometimes is of higher order than considered in this calculation. 



2.2 Pion Lagrangian 

For one-loop calculations of the nucleon magnetic moments, the SU(2) chiral Lagrangian for 
pion fields ir a (a = 1, 2, 3) in the presence of an external electromagnetic field is needed 
only to leading order. With e denoting the unit charge we utilize || 

4™ = JfTr[V^WU + x^U + X tf], (2.1) 

with the chiral tensors 

/ t? 2 i -> _ 

V^U = d^U-iF^U + iUF^ 

X = 2B M + ... (2.2) 



Here f n represents the pion-decay constant (in the chiral limit) and r a are the usual SU(2) 
(Pauli) isospin matrices. As we are only interested in interactions of u and d quarks with 
an electromagnetic background field, the left- and right-handed (axial) vector source terms 
can be identified as 

Fi = F*= e -A,T*. (2.3) 

Furthermore, M. denotes the quark mass matrix for the case of 2 light flavors u, d. The 
connection between the chiral condensate parameter Bo, the non-zero quark masses and the 
resulting non-zero masses for the pion fields will be discussed later. 



2.3 Nucleon and Delta Lagrangians 
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2.3.1 Leading Order Lagrangians 

The leading order Nit and An Lagrangians (diagonal in N and A) required for a calcu- 
lation with explicit pion, nucleon and delta degrees of freedom in the presence of external 
electromagnetic fields are f| ||]: 



(i) 

N 
(1) 



N v [iv ■ D + g A S ■ u] N v 



iv-D lJ 



P ij A 



9i S 



9 \iv Ta 



(2.4) 



where N v corresponds to the non-relativistic spin 1/2 nucleon field and Tjf denotes the non- 
relativistic spin 3/2 delta field with free 4- vector index \i (in Rarita-Schwinger notation) 
and an isovector index i (in isospurion notation) ||. The 4- velocity vector occurs in the 
non-relativistic reduction of the fully Lorentz-invariant chiral Lagrangians, and denotes 
the Pauli-Lubanski spin-vector — details on calculating with non-relativistic chiral effective 
field theories in the so-called "heavy baryon" regime can be found in [12[ . The chiral tensors 
needed for the calculation of the anomalous magnetic moments read 



D. 



i 

2 



u 



with the isospin indices k 



(r M - zV«) 5 ij - ie^Tr (r k rj , (2.5) 
(1,2,3). As we are only working with two light flavors 



u, d, we identify = | as the isoscalar component of the external electromagnetic 



field. The isospin-3/2 projector is defined as usual as 



=3/2 



T"T J . 



As explained 



in section |2.1| , we are explicitly taking into account the finite nucleon-delta mass splitting 
already at leading order, denoted by the parameter A in Eq. ( |2.4j ). Furthermore, qa and g\ 
are the axial nucleon and delta coupling constants. It turns out that g\ is not needed for the 
calculation of the magnetic moments of the nucleon to the order considered in this work. 



2.3.2 Next-to-leading order (NLO) Lagrangian 

The less known NLO nucleon and delta Lagrangians can be found in 0. Here we only 
discuss the terms pertaining to magnetic moments. Defining the chiral field tensors via 



U ] {d,F u R - d u F? - i [F*, F U R ] }VU + VU {d,F u L - d v Fl - i [F%, F U L ] } VU^ 

(2.6) 

we utilize the following representation of the NLO nucleon and delta Lagrangians || 



(2) 
N 



C 



(2) 



N v 



AM 
1 

2M 



[S" t 5"]((1 + K°)f++ 2 + 



N„ 



(2.7) 



[ST, S?\ (D%D k j - DfD*) g, v 



a 6 if+S i i + 2a 7 iV$S i i 



rpu 

3 
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The two couplings n° v and n° s correspond to the (bare) isovector and isoscalar anomalous 
magnetic moments of the nucleon, taken in the chiral limit. Their strength is determined by 
physics which lies outside of the chiral effective field theories. We therefore treat them as 
free parameters to be determined by a fit to lattice simulations of the magnetic moments, 
as discussed in section [|. Likewise, a®, aj are the two corresponding anomalous magnetic 
dipole moments of A (1232). It turns out that these two couplings do not contribute in the 
calculation of the magnetic moments of the nucleon to leading one-loop order — we therefore 
relegate a discussion of these interesting quantities to forthcoming work. 

2.4 Modified NA Transition Lagrangian 

While the leading order ttN and irA interaction Lagrangians discussed in the previous 
section follow the standard rules of chiral power counting and have been used in many 
calculations, we now present a modified version, more appropriate for our purposes, of the 
leading order chiral nucleon-delta transition Lagrangian in the presence of an external 4- 
vector electromagnetic background field A^. The leading order NA transition Lagrangian 
we propose has the form 

4a = Tt[cAwi + <*ifPS*lN v + h.c, (2.8) 

with 

wl = ^Tr (tV) , f£* = ^Tr ( T *f+) . (2.9) 

It involves the axial (transition) coupling ca as well as the (iso-) vector (transition) coupling 
cy, which govern the strengths of the 7riVA and 'yNA vertices in the chiral limit, respectively. 
In section [5] we will discuss numerical estimates for these parameters. It is standard practice 
to include the axial iVA transition in the leading order nucleon-delta transition Lagrangian 
(e.g. see || [13]) via the coupling ca, whereas the leading term of (iso)vector NA transition 
is usually only taken into account at sub-leading order in the iVA transition Lagrangian (see 
e.g. Eq.(112) of Ref.|J). The main reason for this asymmetric treatment of the axial and 
vector NA transition lies in the "standard" counting ( "naive power counting" ) of the chiral 
tensors in powers of the generic mass or momentum scale e: The (pseudo-) vector scales 
as order e, whereas /+ has dimension e 2 . We note that these counting rules on the one hand 
lead to an asymmetric treatment of the vector, axial-vector iVA-transitions by distributing 
them onto the NLO, respectively LO Lagrangians. On the other hand they result in a 
symmetric (NLO) treatment among the magnetic ^NN, '-/NA and 7AA couplings, which, 
for example, is appealing from the viewpoint of the SU(6) quark model. 

In this work we propose the ansatz — displayed in Eq. fl2.8p — to promote the "off-diagonal" 
magnetic NA transition into the leading order NA transition Lagrangian, while leaving the 
corresponding "diagonal" ^NN, 7AA couplings a% in the NLO Lagrangians of Eq.(|2.7p 
as suggested by dimensional analysis*]. In order visualize that the isovector transition struc- 
ture in Eq. (|2.8| ) is part of the leading order Lagrangian, we assign an intrinsic power e" 1 to 

4 One may wonder why we do not propose to promote all three couplings to the respective leading order 



6 



it 



Cy = Cy 1} 



(2.10) 



rendering the structure c v x ' /+ to scale as 0(e), as expected for a leading order Lagrangian. 
In Ref. (U it was argued that this intrinsic scaling of the (dimensionful) coupling cy should 
be made more explicit by introducing a dimension-free coupling cy via Cy~^ = 5y/A. On 
a computational level that prescription is of course equivalent to the structure in Eq.(|2.8|). 
However, given that we are after the quark mass dependence of the magnetic moments, we 
do not want to prejudice as to which dimensionful quantity — be it A or Mn — sets the scale 
in Cy to result in a dimension-free coupling cy. The reason for this caution lies in the fact 
that the two mass scales A and have quite a different intrinsic quark mass dependence 
of their own (see e.g. the lattice simulations given in Ref. [Q), which would then (in higher 
orders of the calculation) seem to alter the chiral extrapolation curve^ depending on which 
parameterization was used for the leading isovector iVA-transition Lagrangian. In order to 
avoid this ambiguity we therefore choose to work with a dimensionful coupling cy = c v , 
keeping in mind that it carries an intrinsic power of e -1 . 

One of the motivations given in Ref. [|J that the magnetic NA transition coupling should 
scale as e _1 , whereas the magnetic NN, A A couplings of Eq. (|2.6|) obey the standard chiral 
counting rules ~ e°, was the well known fact that quark model SU(6) symmetry factors 
relating cy to the isovector anomalous moment of the nucleon typically underestimate the 
7./VA transition strength by as much as 30% |0J. This observation and the fact that 
many observables in pion photo-/electroproduction jl6) as well as in nucleon Compton 
scattering [I7|| crucially depend on a proper treatment of the Ml '-/NA transition in the 
chiral effective field theory provide a physical motivation for a more prominent role of this 
important structure, aside from the formal discussion given above. 

However, the main reason why we insist on having the magnetic iVA-transition to be 
part of the leading order Lagrangian Eq. fl2.8|) does not lie in its mere numerical strength. 
After all the coupling cy — which only represents the leading term of this transition, taken 
in the chiral limit — might be substantially different^ from the physical Ml 7./VA transition 
strength. Even more important from the point of view of the quark mass expansion of 
the magnetic moments is the fact that this operator — as will be discussed in sections 3,f|— 
produces important non-analytic quark mass dependence in the magnetic moments, which 
turns out to be essential for a meaningful chiral extrapolation. Nevertheless one cannot 



Lagrangians. Wc refrain from doing so because for the "diagonal" couplings k®, ae this would lead to the 
peculiar situation that the anomalous ( "Pauli" ) contribution to the nucleon, delta magnetic moments would 
come in at leading order, while the 1 /M suppressed regular ( "Dirac" ) contribution is generated in subleading 
order via the non-relativistic reduction. 

5 The counterterms in the chiral effective field theory do of course not depend on the quark masses. 
However, the problem discussed here arises if one makes the identification between the associated chiral 
limit parameters (e.g. the "bare" nucleon mass M%) and physical quantities (like M N = 0.938 MeV). 
6 The chiral expansion of the Ml 7 TV A transition form factor based on "naive power counting" is analyzed 
. Given that there are three additional counterterms [^] contributing to this transition, we 
do not have good information on the strength of the coupling cy at the moment. 



Refs.pJ, f9 
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proceed at will promoting arbitrary operators into any order of the Lagrangian. One has 
to show explicitly that the resulting effective field theory can still be renormalized. For the 
specific calculation considered here we demonstrate in the course of this paper that this is 
indeed the case. 

To summarize this central paragraph of our work: As a result of our proposal Eq.( p.8|) 
one obtains a symmetric treatment between the axial and the (iso)vector iVA transitions 
and, accordingly, a modified diagrammatic expansion. The consequences of this procedure 
for the case of the chiral expansion of the anomalous magnetic moments of the nucleon 
are the main physics topic of this paperQ. The chiral (i.e. quark-mass) expansion of any 
observable is of course not changed by boosting operators into different orders of the effective 
Lagrangian. However, different expansion schemes can bring in important operators already 
at lower orders in the calculation, thus making a differing scheme more effective. Before we 
can go into the calculation, however, we first have to specify some higher order couplings 
required for a systematic calculations with our modified power-counting. 



2.5 N 2 LO Nucleon Lagrangian 

The standard ChPT calculation of the magnetic moments in the heavy baryon limit gives a 
finite result at leading one-loop order without any counterterm f2|. Allowing for the possi- 
bility that the pion-cloud can also fluctuate around an intermediate spin 3/2 baryon requires 
the introduction of two counterterms proportional to the octet-decuplet mass splitting A to 
be able to renormalize the leading one-loop diagrams in a theory with explicit pion, nucleon 



and delta degrees of freedom based on "naive power counting", as was shown in Ref.pi 
Labeling the two corresponding structures D\ and D 2 , we note that they are related to the 
couplings B 28 , B 2 $ of Ref.|21|| via D\ = B 2 %/ (An fn) 2 , D 2 = B 2 g / (in f n ) 2 . Here we again 



prefer the notation involving dimensionful couplings in order to avoid speculating about the 
underlying mass scale and its inherent quark mass dependence, analogous to our reasoning 
regarding cy in section [2.4| . As explained in Ref. |^TJ the introduction of these two couplings 



is needed for the renormalization of the magnetic jNN vertex function but does not lead to 
observable consequences, as the two structures are quark mass independent. The finite parts 
of these couplings can therefore be utilized to guarantee decoupling of the delta resonance 
in the limit of fixed quark masses and A — > oo for any value of the regularization scale A, 
as will be discussed in section |3|. 

If one now moves on and modifies the "naive power counting" — as we propose in Eq.( [2.8|) — 
one is not surprised to learn that this also leads to consequences in the most generalN 2 LO nu- 
cleon counterterm Lagrangian needed to renormalize the leading-one-loop diagrams. For the 
particular case of our magnetic moment calculations we find — to leading one-loop order — 
that moving the magnetic A A transition into the leading order Lagrangian as proposed in 
Eq.(|2.8|) in general induces four N 3 LO operators with coupling constants D 3 , _D 4 , Ei, E 2 to 



7 Obviously there are interesting applications of this proposal for calculations of electromagnetic scattering 
processes JlO, [L7j in effective field theory which are being explored Q . 
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LEC 




A 


D x 

D 2 
D 3 

D 4 
Ex 
E 2 


*i[S",S*\f+ 

Ai[S»,S»}V$ 

A 2 i[S»,S»]f+ 

i[S»,SP\)#f+ 


-c A c v g A / (27 vr 2 / 2 ) 
+c A c v g A /(3Q7r 2 f 2 ) 



Table 1: N 2 LO counterterms and their /3-functions contributing to the magnetic moments 
of the nucleon to leading one-loop order. 



move down into the N 2 LO Lagrangian^: 

C% ] = N v {D x Ai[S^S v ]f; v + D 2 Ai[S^S v ]V$ 
+D 3 A 2 i [S», 5"] f+ + D 4 A 2 i [S*, S»] V« 

+Ek i [s», s»\ /+ + e 2 i [s», s»\ x « V$ 

\- [[S», ST\ ((1 + 2k*) /+ + 2(1 + 2k° s ) y£>) v°D v + h.c] 



8M 2 



N v , 
(2.11) 



with 



^Tr (v / Z7 t xv / Z7 t + VU X ] ^U) ■ (2.12) 



Note that in the following we work in the isospin limit m u = = rh and therefore only 
need to consider the isoscalar component of x+- Throughout this calculation we utilize 
dimensional regularization and denote the resulting infinities by the quantity L spelled out 
in Appendix All six counterterms then have the generic structure 

Ctf = C i N v O?>N v , (2.13) 

with Cj = C\{X) + Pi 167r 2 L and the associated (3— functions given in Table [TJ. We observe 
that the four counterterms can be separated into a scale (A) dependent finite and an infinite 
part. Once more we want to stress that based on "naive power counting" one would not 
expect to find the local operators associated with D$, D±, E%, E 2 among the terms of the 
Af 2 LO Lagrangian, as their structures x+ an d f~^ v , respectively V$ scale as ~ e 2 each in 
standard counting. In contrast to the coupling cy — which, in the previous section, has been 
attributed an intrinsic power e _1 due to its importance — D3, D4, E\ and E 2 obtain their 

8 In order to construct a complete set of N 2 LO nucleon counterterms required to renormalize all possible 
1-loop graphs involving pions, nucleons and deltas in the presence of arbitrary external fields for the here 



proposed new form of the leading order NA transition Lagrangian Eq.(2.8), one would have to perform a 
new one- loop renormalization analysis for single nucleon processes as for example done in Ref. [p2| , which is 
beyond the scope of this article. 
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intrinsic power of e _1 solely by the requirement that every effective field theory based on the 
most general chiral Lagrangian should be renormalizable, independent of the particular or- 
ganization of the perturbative expansion. Following the notation of section ^]4| we therefore 
write 



£1 

£3 



D 



(-1) . 



E2 
At 



E. 
D 



(-1) 

; 

(-1) 



(2.14) 



rendering the structures 



e[~ 1] X f / 



fJ,V 



and 



E. 



1) (•) y( S ) 



iO nucleon Lagrangian. 



to scale as ~ e , in accordance with the power of the N 

Finally, we comment on the observation that three of the six (3— functions in Table [I] are 
zero, suggesting that three counterterms would be sufficient to renormalize the magnetic 
moments to leading one-loop order. We note, however, that this simplification only occurs 
for our special choice of representation of the chiral field strength tensors given in Eq.( |2.6| ). 
In the (more widely used) conventions of Ref. [12] — which encode the chiral field strength 



tensors as F+ u 



f+ + 2V$ and Tr(F+ ) = AV, 



one would require six counterterms 
with non-zero (3— functions. Another reason for our retaining of D 2 , -D4, E 2 in the N 2 LO 
Lagrangian Eq. (|2.11|) lies in the fact that we want to treat the (unknown) short distance 
physics^ in the isoscalar and the isovector sector of the theory in a symmetric fashion. With 
these remarks we close our discussion on the required chiral Lagrangians and move on to 
the details of the calculation. 



3 The Calculation 

At first glance the reader may wonder why we present yet another calculation for the baryons' 
magnetic moments, as this topic is probably the best studied one in the field of chiral effective 
field theories. The leading non-analytic quark mass dependence is known since the 1970s 
p3[| , calculations with (octet) meson and (octet) baryon degrees of freedom were pioneered 



in baryon ChPT in the late 1980s the advent of heavy baryon techniques brought about 
many more calculations throughout the 1990s, including the first studies with intermediate 
decuplet baryons f2"5 ]. An overview of some recent calculations and references can be found 



111 



When we present our leading-one-loop results for the anomalous magnetic moments 
of the nucleon arising from nucleon, delta and pion degrees of freedom as given by the 
diagrams in Fig.|I] and the corresponding amplitudes in Appendix |A], we do not claim that 
we have calculated any new contributions previously not considered in the literature. In fact, 
Ref.|25|] contains even more one-loop diagrams than we consider here. We point out that the 



underlying philosophy between our work and for example Ref.[^5| is a different one. In ]25 
the authors have calculated all possible one-loop topologies contributing to the magnetic 
moments based on the leading photon and meson couplings to octet/decuplet baryons. In 

9 It also turns out that — for the case of 2 flavors considered here — E2 provides the leading quark mass 
dependence in the isoscalar sector, see section [| 
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this work our power counting — as discussed in section — establishes a hierarchy among the 
one-loop diagrams, selecting the ones given in Eq. ( |A. 1| ) to be leading-one-loop (= 0(e 3 )) and 
dictates the structure of counterterms to be included at this order (c.f. Eq. fl2.1l|) ). Given 
that we spent the first part of this paper arguing for a "modified power counting" it should 
be obvious that chiral effective field theory for low energy baryon properties does not possess 
one unique perturbative expansion parameter like oiqed in Quantum Electrodynamics. Due 
to the complex structure of the low energy hadron spectrum and the many different scales 
involved several counting approaches have to be explored, our proposal in Eq. fl2.8Q is only 
one possibility, albeit a well-motivated one. Ultimately the success in describing phenomena 
determines which expansion scheme really is "effective". With this in mind we point out 
that our leading-one-loop calculation of the anomalous magnetic moments of the nucleon 
depends on only a few structures and parameters, providing some "predictive" power with 
respect to the few lattice data available at the moment. This will be discussed in the 
upcoming subsections. Of course it is possible to push the calculation to one higher order 
(= 0(e 4 )) P7| , which formally includes 35 extra diagrams, containing the remaining ones 
of Ref .12511 as well as some new ones 10 . On the other hand, at 0(e 4 ) one encounters at least 
four additional unknown couplings, making the comparison to lattice data more strenuous, 
so we defer it to a later stage [27| . 



The NLO calculation discussed here proceeds in a straightforward manner, the results 
can be found in Appendix To the order we are working 10 one loop topologies displayed in 
Fig.p] have to be analyzed — of which only 4 yield non-zero results (c.f. Appendix 0). We note 
that the isovector iVA-transition coupling cy discussed in the previous section contributes 
to two of the loop diagrams. We utilize dimensional regularizationF] throughout, the (3- 
functions of the six counterterms are given in Table |I|. The finite parts of the counterterms 
E[(\), -E£(A) will be left as free parameters and fixed from lattice data in section [5]. However, 
as already indicated in section ^75] , matters are different in the case of D\, i = 1, ... 4: The 
infinite parts of these four counterterms are utilized to cancel divergences ~ AL and ~ A 2 L 
(c.f. Table p]). The finite parts of these four structures cannot be observed separately from 
the chiral limit structures k®, k®. As suggested in Ref.pTJ, one can make use of this freedom 



in implementing decoupling of the Delta resonance, e.g we demand that the theory with 
explicit delta degrees of freedom transforms itself into a theory with just pion and nucleon 
degrees of freedom in the limit A — > oo (for fixed quark masses). This implies, for example, 
that all quark mass independent polynomial structures in A must vanish, for any value 
of the chosen regularization scale A. With the results from Appendix A we find that the 
assignment 

D r 2 (X) = 0, 



10 At Oje 4 ) one also has to take into account the first 1/M corrections to the axial NN and NA vertices, 
e.g. see ||. 

n For recent work analyzing the magnetic moments of the octet baryons employing lattice regularization 
methods sec Rcf.[E8f. 
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r>rr\\ c A c v g A Jl , 2A 
Dl(X) = 0, 



(3-1) 



satisfies the condition of vanishing polynomial structures. The decoupling of the Delta 
resonance achieved by this choice can be best seen in the limit m n /A — > 0, which we discuss 
in section |4.2.2|. After these technical comments we finally proceed to the physical results. 



4 Analytic Results 
4.1 General Remarks 

We now turn to the quark mass dependence of the proton and neutron magnetic moments, 
fi p and /i n . Two remarks are in order at this point. First, we note that the chiral corrections 
concern only the anomalous parts of the magnetic moments 

Hp = 1 + K p , jJ n = K n . (4.1) 

Second, the chiral corrections affect the isovector and isoscalar anomalous magnetic moments 
of the nucleon k v , k s quite differently. We therefore discuss our results in the isospin basis 
defined by 

K v = Kp K n , K s = Kp + K n . (4-2) 

Both k v and k s are functions of the isospin averaged quark mass0 rh = (m u + m^)/2 and of 
the chiral condensate parameter Bq. These two quantities are combined to form the leading 
order term in the quark mass expansion of the (squared) pion mass. One finds || 

ml = 2mB {l + O(mB )} 1 (4.3) 

which to leading order corresponds to the well-known Gell-Mann, Oakes, Renner relation 



2Sf| . Changing the value of the light quark mass m (in a numerical simulation of QCD) 
therefore leads to a quadratic change of the pion mass (modulo higher order corrections). 
In the following we present our results as functions of m^. When comparing a magnetic 
moment calculation performed in a chiral effective field theory to lattice QCD data, one 
therefore requires these data as a function of the mass of the lowest lying 0~ boson in the 
simulation, which is identified as the corresponding lattice pion. It is understood, that both 
the pion mass and the associated nucleon magnetic moment simulation are performed with 
identical lattice parameter input. Such correlated lattice results have been reported in the 
literature and we will discuss them in section |5|. Here we focus on the analytic results. 



12 In addition to neglecting all effects from strong isospin breaking — i.e. the m u — md mass difference — 
we also do not consider electromagnetic corrections arising from the different charges of the light u and d 
quarks. All other quarks taken to be infinitely heavy and are effectively integrated out of the theory 
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4.2 Isovector Anomalous Magnetic Moment 



4.2.1 Analytic Results to leading one-loop order 

All results presented in this section can be directly read off from the amplitudes shown 
in Appendix refering to the relevant Feynman diagrams shown in Fig.[I|. . The NLO 
result of SU(2) non-relativistic ("heavy") baryon ChPT with only pion and nucleon degrees 
of freedom (in the following denoted clS CclSG A) reproduces the well-known f2]| leading 
non-analytic quark mass correction to the isovector anomalous magnetic moment of the 
nucleon 

< = <A~ 9 -^ff + N 2 LO. (4.4) 

The standard power-counting of heavy baryon ChPT tells us that to NLO this result orig- 
inates from graphs a) ... e) of Fig.|l} All other contributions are relegated to higher order 
corrections^. 

Recently the chiral corrections to n v have also been evaluated to NLO in the SU(2) Small 
Scale Expansion approach plf , which includes explicit pion, nucleon and delta degrees of 



freedom (in the following denoted as case B). With the additional counting prescription for 
the finite scale A (c.f. section |2.1| ), this calculation follows the standard ("naive") power 
counting rules of baryon ChPT and thus only includes the coupling ca in the NA transition 
Lagrangian of Eq. ( 



B „0 9 A m * M 



B 



2c^AM I L ml 




9^/* V A^ 



2A 



I log [R(m„)] + log }+N 2 LO, (4.5) 



with 




R(m n ) = — + (4.6) 



TTl-rr \l TU 2 



To NLO this result originates from graphs a) ... i) of Fig.|T], all other contributions are again 
relegated to higher orders. 

Now we present the chiral corrections to k v calculated to NLO in our modified scheme 
(in the following denoted as case C). In contrast to the calculation of Ref. [21], we use the 



modified leading order nucleon-delta transition Lagrangian Eq.( |2.8| ) which includes cy — the 
leading term of the isovector 7./VA Ml transition — as well as the induced additional N 2 LO 



13 The same calculation performed to NLO in relativistic baryon ChPT is discussed in Appendix [b]. It 
contains many more structures which are part of the higher order corrections in the heavy baryon approach 
discussed here. 
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counterterm E\ of Eq.( [2.1l| ), required for renormalization of the extra one-loop graphs 
involving cy. We obtain 



K 



C 



g\ m n M 




QP ,u M 2 , 4:C A c v g A Mml 
-8E l {X)Mm n + ^—- 2 log 



log [R(m n )} + log 
2A 



A 



Ac A cvg A Mml 

277T/2A 

8c A c y ^A 2 M 

27tT 2 / 2 



+iV 2 LO . 



ml 



3/2 



2A 



1 - ) log [i^m.)] + 



3m 2 
2A 1 



log 



2AJ 



(4.7) 



This NLO result arises from the graphs displayed in Fig.|l|. We note that very few diagrams, 
of which only 5 are non-zero (c.f. Appendix [A]) to this order, produce such an intricate 
quark mass dependence in k v , making the modified scheme proposed here very effective in 
calculating chiral corrections for electromagnetic quantities. Before we study the quantita- 
tive differences in the chiral corrections to k v between schemes A, B and C, we first draw 
a qualitative picture of the extra physics contained in scheme C. For a proper comparison 
we therefore move on to a discussion concerning the chiral limit. 



4.2.2 Chiral Limit Results 

For completeness we show here the chiral limit results of the three different effective field 
theory calculations for k v discussed above — all calculated to NLO accuracy by Taylor ex- 
pansions of Eqs.( P~4]P~o] , |4T7| ): 



K v \NLO 



K'v \NLO 



K v \NLO 



K v,B 



9\M 



m m 



c\M 



+ml [0 + N 2 LO] + ml 
,n 9aM 



18tT 2 / 2 A 

144vr 2 / 2 A s 



9tt 2 / 2 A ° g V2A/ 
36tt 2 / 2 A 3 & \2AJ 



N 2 LO 



N 2 LO 



+ 



v v,C 



4vr/- 



c\M 



+ 



1 (—\ 

18tt 2 / 2 A ' 97r 2 f 2 A ° S V 2 A / 



8P / UA# . 4c A c v g A M 
—8Ex(X)M + — log 



2A\ 



2c A c v g A M 

277T 2 / 2 



N 2 LO 
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c A cvg A M c A c v g A M /m^N 
12A 2 tt 2 / 2 9A 2 tt 2 / 2 ° g V2A/ 

(4.8) 

First we observe that schemes B and C produce a whole string of terms proportional to 
m™. We note that all these terms — with exception of the chiral limit couplings k® — do 
occur at the same NLO order in the chiral powercounting. They just represent the first 
few terms of an infinite Taylor series in m n arising from the chiral limit expansion of the 
logarithms in Eqs. (|4.5|J4.7| ). We note that all these Taylor coefficients starting from m\ will 
receive corrections at N 2 LO and higher orders. One also observes that the decoupling of the 
Delta resonance as discussed in section |3| is manifest^ in Eq. (|4.8|) , based on the counterterm 
prescription given in Eq.( |3.1| ). 

One can now clearly see that the power counting of scheme C incorporates the first four 
terms in the chiral expansion of k v already at NLO. Specifically, we note that the structures 
proportional to m\ are absent to this ordeiQ in schemes A and B. In the traditional 
expansion scheme A of Heavy Baryon ChPT such terms could only be generated at N 3 LO 
(i.e. at the two- loop level), thus making the explicit calculation of such terms extremely 
prohibitive. The analytic results presented here and in the previous section are completely 
general. In order to produce chiral extrapolation functions for the magnetic moments which 
connect lattice data to the physical worldline and even to the chiral limit we now move on 
into the numerical analysis of the results presented so far. 



i 4 



4:C A c v g A M 
27tt/ 2 A 

144vr 2 / 2 A 3 
+N 2 LO] 



+ N 2 LO 



+ 



c\M /m.N 
36tt 2 / 2 A 3 S V2Ay 



+ 



5 Numerical Results 



5.1 General Remarks 



In principle all couplings and masses — aside from m n — occuring in Eqs. ( |4. 1[|1.5| . [4~7|) are to be 
taken at their values in the chiral limit. However, for most of them the chiral limit values are 
only poorly known. On the other hand, the difference between the values taken at = 
and the physical values of the couplings is of higher order, allowing us to resort to physical 
parameters — given in Table |2] — in those cases where we have only limited information. 

We note that to the order we are working the nucleon mass M is not really a parameter 
occuring in the calculation, as can be easily seen from Eq. ( |A.1|) . Accordingly, an overall 
M can be factored out on the right hand side of Eqs.(4.4,4.5, 4~7|) . The scale M is just a 
convention to obtain magnetic moments in units of nuclear magnetons ([n.m.]) |52]| . For 



14 Thc logarithmic scale dependence ~ log A/ A in Eq.(4 
the delta resonance in the decoupling limit A 
A if A > A x , where A x denotes the chiral symmetry breaking scale. 



retains the information about the existence of 
Decoupling is presumably already achieved for finite 



15 Structures ~ m% are also generated at NLO in relativistic Baryon ChPT, as discussed in Appendix 
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Parameter 


Physical Value 


n a 
iJA 


1.267 


c A 


1.125 


U 


0.0924 GeV 


M 


0.9389 GeV 


A = Re [M A ] - M 


0.2711 GeV 



Table 2: Input parameters used in this work. With "physical value" we denote their magni- 
tudes at the point where the lowest lying Goldstone boson in the theory has the mass of 138 
MeV and is identified with the pion. The physical meaning of these parameters is explained 
in section |2|. The nucleon and pion masses are taken as isospin-averaged. 



the axial coupling constant gA of the nucleon we use its physical value [p2| . Not much 
is known about its value in the chiral limit g% but recent lattice data suggest that the 



quark mass dependence of this quantity is rather flat |33fl . For the pion decay constant we 
utilize its physical value 0, as the difference to /° is known to be be only a few percent 
|§. In order to fix the parameter A — representing the nucleon-delta mass splitting — we 
employ the dispersion theoretical analysis of Ref . [ 33] to obtain the real part of the complex 
delta mass Ma = (1211 — z50) MeV. Recent lattice simulations discussed in |l4j] suggest 
that A also has a rather weak quark mass dependence. Finally, we fix the leading axial 
iVA coupling constant ca by reproducing the imaginary part[^] of the delta mass given in 
Ref. |34J . Unfortunately the quark mass dependence of the axial NA couplings is not known. 
We therefore assume that the physical value for ca constitutes a decent approximation for 
c A . 

If it turns out that some aspects of the reasoning presented here do not hold, then the 
other couplings — which are directly fitted to the lattice data in the upcoming two sections — 
have to compensate for any wrong assignments. We consider this issue to be only a tem- 
porary problem. As soon as more low mass lattice data for a variety of nucleon structure 
operators become available, one can fit all parameters directly to lattice data. 



5.2 Numerical Analysis of Schemes A and B 

We first discuss the numerical results for schemes A and B. With most of the parameters 
of Eqs.( |4.4| , [4.5|) determined from known physical quantities (resulting values are shown in 



Table 0) , it is clear that we have one unknown to fit in each case, n° v A and b . One quickly 
realizes that neither one of the two extrapolation functions provides an m^-dependence that 
is compatible with the lattice data discussed in || . We therefore decide to fit k° v a and k° v b 
in such a way that we reproduce the physical k v = 3.706 [n.m] |32] for m % — > 0.138 GeV. The 
result of scheme A (NLO HBChPT, Eq. (|4.4| )) is shown in Fig.^|by the dashed curve, whereas 
scheme B (NLO SSE, Eq. (|4.5|)) is represented by the dot-dashed curve. While both curves 



16 The value of ca given in Table || corresponds to a strong decay width of Ta = 100 MeV |Q. We note 
that the delta properties given in Ref.[p4| are evaluated at the T-matrix pole. 
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show a rather similar chiral limit behavior, they must be considered inadequate for m 7T > 400 
MeV, as they even change sign in this region^. In fact, the NLO extrapolations shown in 
Fig.Q provide such a poor m^-dependence compared to the lattice data, that a (hypothetical) 
N n LO higher order calculation in scheme A or B which might be more compatible with the 
data shown in Fig.|] would constitute such a large correction to the NLO result presented 
here, that one would have to worry about the "convergence properties" of the respective 
expansion scheme. To summarize our discussion on the numerics for schemes A and B we 
conclude that even future higher order calculations in these two approaches presumably can 
only describe the m^-dependence of the isovector anomalous magnetic moment for m 7T < 400 
MeV, making the comparison with lattice data impossible for the moment. We also note 
that the relativistic version of scheme A discussed in appendix [B| seems to work well out to 
m n « 500 MeV, as shown by the solid line in Fig.|2|. 



5.3 Numerical Analysis of Scheme C 

The situation is more complicated for the NLO result of scheme C (Eq.( [4.7| )), where — at a 
given scale A — we have, in principle, 7 parameters to deal with: 



K 



g° A , /°,c°,A , £[(A),4. (5.1) 



Ideally, we would fit these parameters to lattice QCD simulations of the isovector anomalous 
magnetic moment for < 600 MeVQ. However, at the moment the data situation does 
not allow fits in this parameter range || . Based on our estimate regarding the applicability 
of leading-one-loop calculations in chiral effective field theories we could just stop here with 
our analysis and wait for future simulations at smaller quark masses before we continue 
to discuss the extrapolation curve given by Eq. (|4.7|) . However, we observe that the lattice 
data discussed in Ref.j^j are basically flat or at most weakly dependent on m w in the range 
600 MeV < < 1 GeV~ A x (c.f. Fig.||). This observation leads us to the following 
hypothesis: 

If the extrapolation function of a chiral effective field theory calculation contains sufficient 
quark mass dependent structures to yield such a "plateau" as suggested by the lattice data 
in this mass range, then the higher order corrections — though formally large (~ (m n /A x ) n ) 
form^ > 600 MeV — are presumably small, as they cannot deviate much from the plateau to 
which the extrapolation curve is fitted. 

The "weak" mass dependence of the lattice data for large quark masses in this view acts 
as "boundary condition" for the chiral extrapolation function, constraining the extrapola- 
tion to small quark masses more efficiently than expected from a chiral powercounting for 
individual polynomial structures. Obviously a necessary condition for this hypothesis to 



17 Evcn when loosening our input provided by Table by allowing the (chiral limit) parameters of 
Eqs.(4.4,4.5) to be slightly different from the physical values shown in Table |], the qualitative picture 



in Fig.y does not change. 

18 We take this "naive guess" on the range of applicability of a leading-one-loop calculation based on the 
success of chiral effective field theory for predicting virtual Compton scattering cross sections on a nucleon 
at three momentum transfer \q\ ~ 600 MeV |35[|. 
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make sense is the requirement that the resulting extrapolation curve compares reasonably 
well with the empirical value for k v at m n = 138 MeV. 

With this hypothesis in mind we now leave the secure realm of chiral effective field theory 
and attempt to find a set of numerical values for the parameters in list |5.1| that is consistent 



with the lattice data of Ref.|| in the range 600 MeV < < 1 GeV and still gives a 
meaningful "prediction" for k v at m n = 138 MeV. To pursue this program we employ the 
following philosophy: We determine the number of degrees of freedom that can be fixed from 
lattice data for effective pion masses m n < A x « 1 GeV. For the correlated (k, m n ) lattice 
data given in || it turns out that this number is three[^. We then pick the corresponding 



number of couplings from list |5.1| according to the following principles: 



1) Couplings for which we have no other physical information available. 

2) Couplings for which we expect a significant difference between the physical and the 
chiral limit value. 



Based on the first argument we pick k® c and E{(\) from list |5.1| . This leaves us with one 
more parameter that we can constrain from lattice data. Based on argument two, we choose 
cy, the leading coupling of the magnetic Ml ^yNA transition in the chiral limit introduced 
in Eq.([2.8|). Indeed at the physical point one expects sizable complex- valued corrections to 



this transition from the pion cloud of the nucleon [18, PRl (as well as interference from three 



additional higher order couplings [H5|). The remaining four couplings — g A , c° A , /°, Ao — are 
taken at their physical values gA, ca, f^, A, as given by Table |], i n complete analogy to 
the procedure in schemes A and B discussed in the previous section. The induced quark 
mass difference is again considered to be part of the N n LO, n > 2 corrections in Eq.( |4.7| ). 
Next we utilize the three lattice points of Ref. || together with the parameters given in 
Table ^| as input and generate — at a chosen scaleQ A — numerical estimates for the three 
unknown parameters k® c , E[(X) and cy shown in Table [3]. The resulting couplings are of 
reasonable size^ and produce the chiral extrapolation curve of scheme C shown in Fig.|3|. 
The full curve is obtained by fixing the three couplings from the central values of the lattice 
data 0, whereas the dotted curves indicate the error band resulting from the errors of the 
lattice data [38|]. For large pion masses the curve reproduces — by construction — the nearly 
flat behavior suggested by the lattice data, whereas for low masses (m n < 400 MeV) one 
observes considerable curvature. Surprisingly the curves extrapolate rather well into the 
low mass region, as indicated by the full circle representing the physical k v . 



19 There are actually four lattice points below 1 GeV discussed in ||, based on the simulations of Refs.[p6[ 
|37| . However, the two points around m T = 800 MeV are only separated by 5 MeV, effectively only providing 
one degree of freedom for the fit. 

20 The chiral extrapolation curves shown in FigsJ^-fJ do of course not depend on the choice of A. 

21 It turns out that the radiative decay width of A(1232) estimated from the fitted value of Table || is by 
a factor of 4 smaller than the number given in ||32f . Either this means that the properties of the magnetic 
7-/VA transition are really substantially different in the chiral limit than at the physical point, or that there 
are significant higher order corrections which get lumped into an averaged number for cy . This issue can 
only be decided in an N 2 LO analysis p7]|. 
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Parameter 


A = 0.77 GeV 


A = 1 GeV 




K v,C 

c°v 
E[(\) 


\0 . -L — f— U.t: 

- (2.26 ±0.06) GeV" 1 
-(4.4 ±0.1) GeV^ 3 


4-^ i m n 4 

-(2.26 ±0.06) GeV" 1 
-(3.85 ±0.1) GeV" 3 


K s,C 

E r 2 {\) 


-0.11 

±0.074 GeV" 3 


-0.11 

±0.074 GeV" 3 



Table 3: Values of the three isovector and two isoscalar parameters obtained from fitting to 
the lattice data of Ref.0 for different values of the regularization scale A, using as additional 
input the parameters displayed in Table |2|. 



One can now ask the question "how large is the isovector anomalous magnetic moment 
k v in a world where the lattice pion mass is 138 MeV ?" Given that this piece of information 
was not used in determining the parameters in Tables |2| and || we obtain — via Eq. fl4.7|) — the 
"prediction" 

«i;|m„->138MeV = 3.5 ± 0.4 [n.m.] . (5.2) 

This does not seem to be a great achievement, given the experimental accuracy to which 
the anomalous magnetic moments of proton and neutron are known (k% cp ' = 3.706 . . . [n.m.] 



32]|). However, it is by no means obvious that our chiral extrapolation should come any- 
where close to the experimental number in view of the rather large extrapolation range, the 
sizable error bars of the lattice data, the non-negligeable curvature required for a successful 
extrapolation to small quark masses and the large associated mass scales of the lattice pion. 
We therefore consider the result of Eq.( |5.2| ) a rather surprising success. Along the same 
lines one can also determine the isovector anomalous magnetic moment of the nucleon in 
the chiral limit. We obtain the prediction 

«v | tti^-^o = 5.1 ± 0.4 [n.m.] , (5.3) 

i.e. we find a significant enhancement over the value at the physical point. We note that it 
is the coupling of the nucleon's pion cloud to the external electromagnetic field applied to 
probe the strength of the nucleon's magnetic moment, that is responsible for this reduction 
in k v when one slowly increases the quark masses from to about 8 MeV. Diagrammatically 
this effect is best displayed by diagram d) in Fig.jl], which provides the bulk of the effect 



and corresponds to the Caldi-Pagels term displayed in Eq.(4.4). 

Ultimately the hypothesis formulated above about the applicability of Eq. (|4.7|) into the 
realm of m n > 600 MeV due to suspected strong cancellations among the higher order 
corrections can only be tested once the N 2 LO corrections are fully calculated and analyzed 
P7| , but the non-trivial physical and chiral limit predictions given by Eqs.( [5.2| , |5.3| ) look 
rather promising in this respect. We now proceed to the (numerical) chiral limit discussion, 
which will provide some insight into the dynamical origin of the successful extrapolation 
function of scheme C generated in this section. 
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5.4 Numerical Analysis of the Chiral Limit 



In this section we discuss the interplay between the various contributions to k v by looking 
at the chiral limit expansion of all three schemes discussed in section |4.2.2j . With Tables 
and 131 we obtain 



K^[n.m.}\ NLO = 5.645 - -Q^T m * 



14.05 



14 05 TTlP 1 / 771 \ 

*[n.m.]\ NLO « 5.859 - — m?r + ^ (0.6480 + 5.780 log 



k v [n.m.\\ NLO 



4 

mz 



^(16.95 + 19.66 log ^)+... 



GeV 4 V °GeV, 

14 05 717^ / 771 \ 

K°[n.m.]\ NLO « 5.109 -■^ 7 m» + ^i J (36.63 + 5.780 log 

- 61 - 63 ^ + ^( 2443 - 34 61l0 ^) + -' (54) 
which leads to the following observations: 

1) These approximate formulae only hold for pion masses below 400 MeV. For detailed 
numerical studies one should use the formulae of Eq. ( |4.7| ) which contain the full ana- 
lytic structure. 

2) These formulae display Taylor coefficients of a series in m™. There is no hierarchy in 
these numbers in the sense that coefficients for small values of n are larger than the 



ones for high values of n. As already explained in section 4.2.2, these coefficients just 



arise from Taylor expanding the logarithms of Eqs.( |4. j| . [1.5| . [4.7|) . All the structures- 



except for the quark mass independent leading terms — displayed in Eq. (|5.4|) are in 
fact part of the same chiral order (NLO) in their respective expansion schemes A, B 
and C. It is therefore not meaningful to look for "convergence" in this representation 
of the quark mass expansion of the magnetic moments. Convergence will be studied by 
calculating the N 2 LO corrections to Eqs. ([4.4|J4.5| , [1.7|) , Taylor expanding the full result 
as done here and then comparing by what amount the individual Taylor coefficients 
displayed in Eq. (\5.4\) have changed when one moves from NLO to JSFLO accuracy. At 



present the N 2 LO analyses in schemes B and C do not exist yet [27]. As indicated in 
section [5] they involve several new unknown couplings, therefore a detailed evaluation 
can only take place if more lattice data (preferably at smaller quark masses) become 
available in the correlated (k, m^) representation employed here. 



3) Eq.( |5.4| ) gives a good idea why scheme C is so much more effective than scheme B 
with its "naive power counting" . The extra structure ~ m\ actually carries a rather 
large coefficient^ and also the terms ~ m\ are significantly enhanced in C. As can 
be seen from Eq. (|4.8|), both structures are intimately connected with the isovector 



22 For a discussion of these structures in leading-onc-loop relativistic baryon ChPT we refer to Appendix 

I 
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N A transition governed by cy and the induced counterterm E\ . Given that these two 
couplings have such a strong impact on the quadratic and the cubic term in the chiral 
expansion, our proposal of Eq. fl2.8Q seems well justified. 

With these remarks we close our discussion on the chiral extrapolation function of the 
isovector anomalous magnetic moment and move on to the isoscalar sector. 

6 Isoscalar Anomalous Magnetic Moment 

In contrast to the isovector anomalous magnetic moment, there is hardly any quark mass 
dependence for the isoscalar anomalous magnetic moment of the nucleon at NLO in chiral 
effective field theory calculations. In the three schemes discussed above one finds 



Only scheme C shows any dependence on the quark masses to this order |4]]. We observe 
that this quark mass dependence is not related to chiral dynamics but solely arises from 
an internal quark mass dependence of the core/bare spin 1/2 nucleon. Its origin (and the 
strength of the associated counterterm E 2 ) is therefore outside the range of the effective 
field theory. 

An extrapolation of the isoscalar moment data to smaller quark masses analogous to 
the one in the isovector sector is therefore unreliable at this order. Moreover, due to the 
smallness of the isoscalar anomalous magnetic moment and the correspondingly large error 
bars of the lattice simulation reported in [[3] , we must conclude that at present extrapolations 
in the isoscalar sector are not feasible. Fig.|] summarizes the lattice data situation for pion 
masses below 1 GeV. The curve shown in this figure is obtained by fitting k° sC and E 2 to the 
physical value and to the lattice point at m n = 800 MeV. The resulting fit parameters are 
given in Table |3|. Our result in the isoscalar sector therefore is not meant as an extrapolation 
but at the moment merely serves as a "best guess" (within scheme C) for the quark mass 
dependence in k s . Nevertheless we have to utilize this result later when we compare with 
the Pade fits of the Adelaide group for the magnetic moments of proton and neutron. 

Finally, we briefly comment on the leading quark mass dependence to the isoscalar 
anomalous magnetic moment of the nucleon, that does arise from chiral dynamics. It is 
known that at the one-loop level one obtains a contribution if one extends the theory to 
the case of three active flavors — SU(3) Baryon ChPT — with quark masses (fh,m s ). For the 
leading non-analytic quark mass dependence in a theory with only octet mesons and octet 
baryons as active degrees of freedom — corresponding to scheme A in the SU(2) sector — one 



A 



k° S)A + + N 2 LO 
nl B + + N 2 LO 
K° s C -%E 2 Mml + N 2 LO . 



s 




(6.1) 



finds |3 



SU{3) 




(5D 2 - 6DF + 9F 2 ) + N 2 LO 
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Proton 


Value 


Neutron 


Value 




3.31 [n.m.] 


u° 

H'n 


-2.39 [n.m.] 


Olp 


1.37 GeV- 1 


Ql n 


1.85 GeV" 1 




0.452 GeV^ 2 


Pn 


0.271 GeV" 2 



Table 4: Parameters used in the Pade fit of the magnetic moments, taken from Ref. ||. 
with 

m \ = (m + m s )B {l + 0(M)} (6.3) 

denoting the kaon mass (squared). We note that this contribution arises from diagram d) 
in Fig.|l] if one allows the complete baryon and meson octet as possible intermediate states. 
F n = (fn + /if)/2 denotes the average of the pion and kaon decay constants to this order and 
F, D are the SU(3) axial coupling constants |§ with the constraint g\ = F + D. We also 
note that for the case of three active flavors the isoscalar anomalous magnetic moment of 
the nucleon arises from two short distance counterterms {e.g. see the Lagrangian discussed 
in fl2"5|). Given that the present lattice data in the isoscalar sector as shown in Fig.[| are not 
sufficiently accurate, we do not embark on a numerical study of Eq.( |6.2| ) or its combination 
with Eq.Q. 



7 Comparison to Pade Approximants and the Quark 
Model 

In 1998 the Adelaide group suggested || a simple parameterization for the quark mass 
dependence of the nucleons' magnetic moments based on Pade approximants. For their 
"best fit"Q they used the following functional dependence, 

_ Pp,n /rr i \ 

f^p,n — .J n 2 ' I'- 1 / 

with values of the parameters given in Table |4]. In Fig.[5] we show this Pade fit as the dashed 
curve, whereas our combined central value result in scheme C 



if v = 1 + {k c s +k^)/2 



p 

Mn = K C -0/2, (7-2) 

is represented by the full curve, utilizing the parameters given in Tables ^| and ^ for k% 
from Eq.(4/7) and from Eq.( |6.1|) . Surprisingly, both parameterizations agree quite well 



and are certainly compatible within the present error band. Minor differences can be found 
in the curve for the neutron and for really small pion masses near the chiral limit, but 



23 The Adelaide group also includes lattice points above 1 GeV pion mass in their analysis and obtained 
a good fit throughout the whole region in m v . 
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at present both parameterizations are indistinguishable due to the sizable extrapolation 



errors originating from the lattice error bars p8| . Both parameterizations of the chiral 
corrections to the magnetic moments indicate that for small quark masses approaching the 
"physical" values, there is indeed a substantial curvature with respect to beyond the 
leading-non-analytic Caldi-Pagels term of Eq. (|4.4j ). We note again that, a priori, there is 
no reason to expect that a simple Pade ansatz as in Eq. (|7.1|) captures all the important 
chiral physics, especially in view of potentially large logarithmic terms generated by the 
pion-cloud dynamics. Judging from the chiral limit result of our chiral extrapolation C 
displayed in Eq.( |5.4| ) the numerical coefficients in front of the m^logrn^ structures are 
either small (c.f. n=l) or get canceled by corresponding m 2 T n polynomial structures (c.f. 
n=2). The dynamical origin of this logarithmic suppression is not known at the moment. 
Nevertheless, we observe that our microscopic calculation agrees well with the extrapolation 
formula Eq. ( |7. 1| ) for the proton and neutron magnetic moments. 

It is instructive to compare the quark mass dependence of the ratio of magnetic moments 
with SU(6) quark model predictions. In Ref. |3{| it was noted that for a pion mass of 
~ 240 MeV the proton to neutron ratio \i v j ' \x n would yield the quark model prediction of 
-3/2, leading to the conclusion that the good agreement between the quark model — which 
knows nothing of the light current quark masses discussed here — and the experimental ratio 
l^p/Hn — —1.46 is accidental. In Fig|| we show that our parameterization for the ratio of 
the magnetic moments given via Eq. ( |7.2p (full curve) follows the trend of the Pade formula 
(dashed curved) Our best fit curve already gives a ratio of -3/2 for m n ~ 150 MeV, 

albeit the error band of our extrapolation (dotted curves) arising from the lattice data is 
sizable in this quark mass regime, pointing again to the need for lattice data at smaller 
quark masses to set stronger constraints on the chiral extrapolation. The restoration of the 



quark model results in the limit of very heavy quark masses discussed in [39 is beyond the 



realm of applicability of chiral effective field theory, so we will not discuss it here. 

While the fi p / fi n ratio shown in Fig.^| allows extrapolation curves to lie near the SU(6) 
quark model prediction at least for certain values of light quark masses, this is not the case 
for the ratio of the anomalous magnetic moments shown in Fig.|7|. The SU(6) quark model 
value K P I K n = —1 is not reached by any of the two extrapolation curves in the small quark 
mass regime. However, it is interesting to note that this ratio seems to be rather insensitive 
to quark mass effects beyond the leading-non-analytic (Caldi-Pagels) term given in Eq.( [4.4| ) 
for pion masses up to 250 MeV, indicated by the dot-dashed curve in Fig.[7] corresponding to 
scheme A of Eqs. fl4.4j , |6.ip . Whereas in the other observables (e.g. see Figs.|3],[5]) the leading- 



non-analytic term is not even sufficient to extrapolate from the chiral limit to the physical 
pion mass, the ratio k p / n n seems well suited for future chiral extrapolations of lattice data 
with small quark masses. With this observation we move on to a brief discussion regarding 
the influence of quenching. . 
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8 Effects of Quenching 



The lattice data on magnetic moments shown in Figs.|2|-|7| have been obtained in the quenched 
approximation of QCD, i.e. with the (loop-) effects of sea quarks effectively suppressed. A 
fully consistent treatment of the chiral extrapolation should therefore be adapted to this 
situation, in the sense that the effects of quenching should also be taken into account 
appropriately within a chiral effective field theory. Such a framework has been developed 
in the past decade, called Quenched Chiral Perturbation Theory (QChPT) pUf . When 
performing chiral extrapolations of quenched lattice data, one does not only have to take into 
account that the prefactors of pion-mass dependent terms could be different in QChPT — 
even the chiral singularity structure can be different between "Quenched" and "Full" QCD. 
Specifically, for the case of the nucleon magnetic moments of interest here, this means 
that there are terms proportional to const, x logm^ which do not exist in "Full" QCD. 
The chiral expansion of k v of Eq.( |4.8| ) demonstrates that the leading logarithmic pion- 
mass dependence is proportional to x logm^ in "Full" QCD. The presence of a new 
chiral singularity like const, x logrn^- in "Quenched" QCD indicates that its low energy 
properties — not to speak of its chiral limit — can differ qualitatively from what we know in 
hadron phenomenology. 

While "Quenched" QCD is interesting in its own right, with the comments just made 
one might arrive at the pessimistic conclusion that not much can be learned about hadron 
properties in "Full" QCD via extrapolations of quenched lattice data. However, the loop- 
effects of the sea-quarks — which are missing in "Quenched" QCD — get strongly suppressed 
for large quark masses, leading one to expect that "Quenched" and "Full" QCD do not differ 
much in the "heavy quark regime". In addition, large quark masses mean that one is "far 
away" from chiral singularities which dominate the chiral limit behavior of the quantities of 
interest. For the time being we can therefore identify quenched lattice data approximately 
with "Full" QCD for effective pion masses above 600 MeV. This assumption is further 
supported by three observations: 



1) Available quenched lattice data on the magnetic moments of the nucleon for pion 
masses in the range 0.6 GeV < m n < 1.5 GeV show a rather moderate curvature (e.g. 
see Fig. 5 in Ref. [§]), suggesting a negligeable effect of the extra chiral singularity in 
this mass region. 



2) The difference between the Pade extrapolation curve Eq. (|7. 1|) discussed in section [7] 
and a similar formula that explicitly includes the leading effects of quenching has been 
reported to be small [^2] . For pion masses above 600 MeV the two analyses lie within 
the error bars of the lattice data. 



3) No significant differences between quenched and fully dynamical simulations for m n > 
600 MeV have been reported for a variety of nucleon structure properties, see for exam- 
ple the recent study of moments of nucleon quark distributions by the LHPC/SESAM 
collaboration [43| . 
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In essence, while "Quenched" QCD does have a chiral limit different from "Full" QCD 
we can accept the given lattice data points for proton and neutron magnetic moments at 
m n > 600 MeV as if they were "unquenched" , within their uncertainties. An extrapolation 
using standard (rather than quenched) effective chiral field theory seems therefore justified 
to provide the correct chiral extrapolation to small quark masses as well as the proper chiral 
limit, at least for our present purpose. 



9 Conclusion and Outlook 

The present analysis has pointed out the feasibility of systematic chiral extrapolations of 
nucleon magnetic moments from lattice QCD, down to the range of realistic light quark 
masses where comparisons with the actual observables can be made. An important element 
in this discussion is the treatment of the A(1232) isobar as an explicit degree of freedom 
in view of its important role in the magnetic structure of the nucleon. In our approach it 
is this feature which produces the important non-analytic quark mass dependencies of the 
magnetic moments beyond the well-known Caldi-Pagels term proportional to m n . 

Our resulting extrapolation is remarkably close to the Pade approximant parameteri- 
zation of the Adelaide group. A sign of caution should be added, however. The existing 
(quenched) data terminate around an effective pion mass of 0.65 GeV, corresponding to 
u- and c?-quark masses rn u d ~ 200 MeV and larger. Expanding one-loop chiral effective 
field theory to such large quark mass scales has its inherent uncertainties which induce a 
substantial error in the extrapolation down to small quark masses. In fact, it is rather sur- 
prising that our NLO extrapolation curve ends up near the physical value for k v , given that 
we only use the input from these rather large mass scales. Certainly the N 2 LO corrections 
in schemes A, B and C (as well as in the relativistic approach discussed in Appendix ||) 
need to be analyzed systematically to judge the stability of the extrapolation fl2~7| . Future 
(partially) unquenched lattice simulations aiming for effective pion masses around 300 MeV 
are also expected to reduce the extrapolation uncertainties significantly. According to our 
results deviations from the nearly linear trend seen so far in the data below 1 GeV pion 
mass should then become visible in the lattice data. However, extrapolating the results of 
such intermediate mass scale simulations down to the "physical" pion mass of 138 MeV will 
presumably require precision calculations within (partially) quenched chiral effective field 
theories, in order to get control over the magnitude of the effects of (partial) quenching in 
the simulations. Pioneering studies in this direction have already been performed and look 
promising 



1J 



As a final remark we note that future lattice simulations (and the associated chiral 
extrapolations) of magnetic moments and related nucleon properties should preferentially 
be done in the isovector/isoscalar basis, as the two channels show quite different patterns 
of quark mass dependence. 
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A Amplitudes 

Here we present the results for the 11 leading-one-loop (i.e. 0(e 3 ) diagrams shown in 
Fig.Q which can contribute to the anomalous magnetic moment of a nucleon of mass M^. 
The Lagrangians needed for this calculation are discussed in section |2|. We work in the 
Breit-frame and choose the velocity vector = (1,0,0,0). With denoting the Pauli- 
Lubanski spin-vector and & denoting the polarization 4-vector of an incoming photon with 
4-momentum one finds 
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Explicit expressions for the function J 2 are given in [|T^|. We evaluate the amplitudes in 
(^-dimensions with induced regularization scale A. Any ultraviolet divergences appearing in 
the limit d — > 4 are subsumed in 



L 



167T 2 



d 



_ + -( Tfi -l + l n 47r) 



(A.2) 



where 7^ denotes the Euler-Mascharoni constant. 

To simplify the calculation we have utilized the electromagnetic gauge-condition v ■ e = 
0. Amplitudes e) and i) are zero due to this choice of gauge, whereas the null result of 
amplitudes b), c), f) and g) follows from the Pauli-Lubanski condition S ■ v = 0. To this 
order in the calculation the non-zero results therefore arise solely from amplitudes a), d), 
h), j) and k). However, diagrams b), c), e), f), g) i) will start contributing^] at C(e 4 ) with 
a strength depending on our choice of u M . 



B 



to NLO in Relativistic Baryon ChPT 



Expansion schemes A, B and C discussed in the main text are non-relativistic approaches 
based on the "Heavy Baryon" method of ||. Recently, the relativistic one loop analysis of 
the nucleon magnetic moments presented in [24] has been updated in Ref.|3(|, employing 
a new regularization scheme JD]] which overcomes the large renormalization effects ~ M n 
typically plaguing relativistic Baryon ChPT in standard dimensional regularization p4 |. At 
leading-one-loop order (i.e. NLO) in a relativistic chiral effective field theory with pions and 
nucleons as the degrees of freedom one obtains [|30] 



JR 



g 2 A M 2 x 2 



c 6 



2 



2 (8 - 13x 2 + 3x 4 1 



arccos 



or 



2 (7 - 3x 2 ) log 



+N 2 LO 



(b.i; 



with x = m^/M. This result arises from diagrams a) -e) in Fig.0, as well as tadpole and 
wavefunction renormalization. We are interested in this result because relativistic correc- 
tions to the Caldi-Pagels term given in Eq.( |4.4| ) can also generate structures ~ m^, which 
are claimed to be important for the fitting of the lattice data (c.f. section |5.4j) . To discuss 
this structure we perform the chiral limit expansion and obtain 



K 



IR\ 

v \NLO 



C6 



9aM 



^ 2 f 2 



+ M! lo ^ + N2LO 



8nMf 2 



+ N 2 LO 



(B.2) 



24 Note thatatO(e 4 ) one also has to take into account tadpole topologies and wavefunction renormalization 
graphs, which are not displayed in Fig^. 
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which agrees with Ref.[24], where terms up to m^logm,,- were discussed. Indeed one also 
finds a structure ~ m\ in the relativistic approach, as expected. In order to judge its im- 
portance we utilize the parameters of Table H and fix the unknown constant c 6 to reproduce 
k v = 3.706 [n.m.] for m n — > 138 MeV. We obtain the solid curve shown in Fig.f2|. Its 
m^-dependence is superior to scheme A or B when compared to the lattice data of Ref. 0. 
However, the NLO relativistic result of Eq.( B.l ) breaks down for m n ps 600 MeV. Returning 



to the chiral limit discussion we find 

14.05 ml 



ZL Mr) 777 / 777 \ 

n.m]\ NLO « 5 - 068 --^ m - + ^(- 5 - 814 - 16 - 671 °g^) 



+23.90- " 



Comparing this result to Eq.( |5.4| ) we conclude that the relativistic approach like scheme 
C contains all allowed quark mass structures, albeit with different (in general smaller) 
coefficients. We note that the structures ~ m^, m\ have a different sign in the relativistic 
approach than in Eq.( |5.4j) . At N 2 LO in relativistic Baryon ChPT some effects related to 
A(1232) are implicitly incorporated^ in the higher order counterterms of the delta free 
theory. It will be interesting to see whether at that order the relativistic approach can be 
extended to pion masses above 600 MeV. A N 2 LO stability analysis of schemes A, B, C and 
of the relativistic approach is in preparation [^7] . 
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Figure 1: Diagrams contributing to the anomalous magnetic moment of the nucleon at 
leading-one-loop order. 
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Figure 2: Pion mass dependence of the isovector anomalous magnetic moment in nuclear 
magnetons. The curves shown denote the standard NLO Heavy Baryon ChPT (dashed line, 
scheme A, see Eq.( [4.4|) ) and the NLO Small Scale Expansion calculation (dot-dashed line, 
scheme B, see Eq. (|4.5| )). The solid line denotes the LO Relativistic BChPT result discussed 
in Appendix [FJ. The lattice data are taken from Ref.[[|. The physical k v = 3.706 [n.m] is 
displayed by the full circle. 
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Figure 3: Pion mass dependence of the isovector anomalous magnetic moment in nuclear 
magnetons. The full curve denotes the NLO calculation in the modified expansion scheme 
C of Eq.( [4.7| ). The lattice data are taken from Ref.||. The physical n v = 3.706 [n.m] is 
displayed by the full circle. 
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Figure 4: Pion mass dependence of the isoscalar anomalous magnetic moment in nuclear 
magnetons. The full curve represents the suggested pion mass dependence of the modified 
expansion scheme C of Eq. fl6.ip . The lattice data are taken from Ref. ||. The physical 
k s = —0.1202 [n.m] is displayed by the full circle. 
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Figure 5: Pion mass dependence of the magnetic moments of proton (upper curves) and 
neutron (lower curves) in nuclear magnetons. The full curve represents the best fit in 
the modified expansion scheme, whereas the dashed curve denotes the Pade extrapolation 
formula Eq. (|7.1|) . Our error estimate is given by the dotted curves. The lattice data are 
taken from Ref. ||. The physical values fi p = 2.793 [n.m], \i n = —1.913 [n.m] are displayed 
by the full circles. 
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Figure 6: Pion mass dependence of the ratio of the magnetic moments. The physical value 
Vp/ fi n — —1-46 is denoted by the full circle, whereas the SU(6) quark model prediction 
of -3/2 is given by the horizontal straight line. The chiral extrapolation result based on 
Eq.( |7.2| ) is given by the solid curve, with extrapolation errors indicated by the dotted lines. 
The corresponding curve based on the Pade formula of Eq. ( [7.1| ) is shown as the dashed line. 
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Figure 7: Pion mass dependence of the ratio of the anomalous magnetic moments of proton 
and neutron. The SU(6) quark model prediction of -1 is given by the straight horizontal 
line. The chiral extrapolation curve based on scheme C is given by the solid and the Pade 
fit by the dashed curve. The physical value n p /n n = —0.94 is denoted by the full circle. 
The dot-dashed curve shows the dependence of this ratio on the leading-non-analytic quark 
mass term only, corresponding to scheme A discussed in the text. 
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